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Within the context of a simple mechanical model the paper 
examines the movement of a satellite with respect to the cente 
of masses under conditions of uniaxial aerodynamic attitude 
control. The equations of motion of the satellite take accour 
of the gravitational and restorative aerodyn5unic moments. It 
is presumed that the aerodynamic moment is much larger than * 
the gravitational , and the» motion eguations contain a large • 
parameter. A two-parameter integrated surface of these 
equations is constructed in the form of foinmal series in 
terms of negative powers of the large parameter, describing 
the oscillations and rotations of the satellite about its 
lengthwise axis, approximately oriented along the orbital 
tangent. It is proposed to treat such movements as 
nominal undisturbed motions of the satellite under conditions 
of aerodynamic attitude control. A numerical investigation 
is made for the above integrated surface. 
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Annotation 


Within the context of a simple mechanical model the paper 
examines the movement of a satellite with respect to the center 
of masses under conditions of uniaxial aerodynamic attitude 
control. The equations of motion of the satellite take account 
of the gravitational and restorative aerodynamic moments. It 
is presumed that the aerodynamic moment is much larger than the 
gravitational, and the motion equations contain a large para- 
meter. A two-parameter integrated surface of these equations 
is constructed in the form of formal series in terms of negative 
powers of the large parameter, describing the oscillations and 
rotations of the satellite about its lengthwise axis, approxi- 
mately oriented along the orbital tangent. It is proposed to 
treat such movements as nominal undisturbed motions of the 
satellite under conditions of aerodyneunic attitude control. 

A numerical investigation is made for the above integrated 
surface. 
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UNIAXIAL AERODYNAMIC ATTITUDE CONTROL 
OF ARTIFICIAL SATELLITES 

V. V. Sazonov 

The Keldysh Institute of Applied Mathematics, 
USSR Academy of Sciences 


1 . Introduction 


The aerodynamic moment can be used on circular and slightly /3 * 
elliptical orbits in the altitude range of as much as 500 km 
for attitude control of the longitudinal axis of an artificial 
earth satellite in regard to the oncoming air stream, the 
velocity of which is practically directed along the orbital 
tangent. For this purpose, the shape of the exterior shell of 
the satellite is chosen so that, when the normal attitude con- 
trol is disturbed, pitch and yaw restorative moments are created 
that strive to line up the longitudinal satellite axis with the 
velocity vector of the oncoming current. Since the projection 
of the aerodynamic moment on the longitudinal satellite axis 
is equal to zero,^ in itself the aerodynamic moment can only 
provide a unieocial attitude control of the satellite. 

This attitude control principle was used aboard the Soviet 
satellites Kosmos-149 and Kosroos-320 [2-4]. However the atti- 
tude control system of these satellites was not a purely aero- 


^For certain configurations of the outer satellite shell the 
projection of the aerodyncunic moment on the longitudinal axis 
is different from zero. This projection may result in a 
twisting of the satellite (propeller effect [1]) and somer- 
saulting. Therefore the shape of the outer shell of a satellite 
that functions under conditions of aerodynamic attitude control 
should be chosen so that this projection is identically equal to 
zero. 

* 

Numbers in the margin indicate pagination in the foreign text. 


1 


dynamic one. In order to create a restorative moment about the 
longitudinal axis, two-degree gyroscopes were Installed aboard 
these satellites, also serving as absorbers of the perturbing 
motion. By using the gyroscopes and the gravitational moment 
It was possible to obtain a trlaxlal attitude control of these 
satellites. But If It Is enough for the satellite to have a 
uniaxial attitude control In terms of the orbital temgent to 
accomplish Its flight mission, the principle of aerodynamic 
attitude control In the pure form can be used. The present work 
Investigates the possibilities of such a method of attitude 
control. 

The satellite discussed below Is regarded as a solid with 
outer shell In the form of a sphere. The center of the sphere 
does not coincide with the center of masses of the satellite 
and lies on one of Its main central coces of Inertia. Such a 
satellite may serve as a model of that shown In Fig. 1. Here 
the Instrument package has small geometrical dimensions, emd 
we can Ignore the Influence of the atmosphere on It. The sphere ^ 
acts as an aerodynamic stabilizer and has a mass much less than 
that of the overall satellite. 

The equations of motion of the satellite take account 
of the gravitational and restorative aerodynamic moments. It 
Is presumed that the aerodynamic moment Is much larger than the 
gravitational, and the motion equations contain a large para- 
meter. A two-parameter Integrated surface of these equations 
Is constructed In the form of formal series In terms of negative 
powers of the large parameter, describing the oscillations and 
rotations of the satellite about Its lengthwise axis, approxi- 
mately oriented along the orbital tangent. It Is proposed to 
treat such movements as nominal undisturbed motions of the 
satellite under conditions of aerodyneunlc attitude control. A 
numberlcal Investigation Is made for the above Integrated sur- 
face. 
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Let us consider a satellite (solid ) , the center of masses 
of which is moving in a circular orbit about the Earth. To 
write the equations of motion of the satellite relative to the 
center of masses we shall introduce two clockwise Cartesian 
systems of coordinates. 

Ox^,X2»x^ - a coordinate system rigidly bound to the 
satellite. The point 0 is the center of masses of the satellite. 
The axes Ox^fOx^^Ox^ are the major central axes of Inertia of 
the satellite. 


0X^X2X^ - the orbital system of coordinates. The axis 
OX^ is directed along the radius vector of the point 0 with 
respect to the center of Earth, the axis OX^ is directed along 
the orbital tangent in the direction of movement of the satellite. 


We shall assign the position of the coordinate system 
0x^X2X^ relative to the system 0X^X2X^ by using the angles 
a, $, Y (Fig. 2 ). The matrix for conversion from one of these 
systems to the other has the appearance: 


/ 

Xi 


*» 


«« 





aii 



a,i 


a» 


^ Und UnX, 


(I ij 

a j 2 - V. 

a,2 «: ayisinli. 


S 

iti\ s aindccsp, 

o $ihdSihj\C4*sX 

(1^3 Sind sinfiUnlf’f^co$dC4fsif. 


We introduce the definitions; ,u)2 # 0 )^) " the ed^solute 

angular velocity of the satellite (hereafter the components of 
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the vectors are given in the coordinate system Ox^X 2 X^) ; ^ 

A,B,C -* the moments of inertia of the satellite with respect to 
the axes Ox^,Ox 2 »Ox 2 ; (M^fM 2 /M 2 ) *■ the major moment of the 
external forces acting on the satellite; Uq « const > 0 <- is the 
angular velocity of orbital motion; t is the time. The movement 
of the satellite relative to the center of masses is described 
by equations: 


(C - - Ml. 

Mi. 

Ms. (1) 

uf * freest. 

cosy 

mU)^ 4 CCS<i * Mi sinnj - Me 7 ^ • 

^ = -ou^ sinti. 4 Wifcs>A * Me sin K. 


Of the external moments acting on the satellite we shall 
only consider the gravitational and restorative aerodynamic 
moments. The components of the gravitational moment have the 
form: 

Mji ' 3u}f(C-B)iluaji , Mjt • 3u)^(A- C)o i^tt a, = . 


In calculating the restorative aerodynamic moment we shall 
consider that the outer shell of the satellite has the shape of 
a sphere of radius R with center at point: 

(«/ , 0 ). 


Here d < 0. With respect to the interaction between the 
satellite and the atmosphere we assume the following: 1) the 

atmosphere is immobile in absolute space; 2) the action of the 
atmosphere on the satellite reduces to a resistance force applied 
at the center of pressure and directed opposite the velocity of 
the center of masses of the satellite; 3) the resistance of the 
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atmosphere does not affect the evolution of the orbit. Under 
these assumptions the components of the aerodynamic moment are: 


Here Is the coefficient of aerodynamic resistance, p Is the 
density of the oncoming air stream, V Is the velocity of the 
center of masses of the satellite. 


Until now the orbit of the center of masses of the satellite 
was taken as circular. Such an assumption Is perfectly justified 
when the eccentricity Is <0.01, although the calculation of the 
density of the oncoming air stream should allow for even a small 
orbital elllptlclty. We shall assiame that this density depends 
on the height of the satellite h above the surface of the Earth 
as : 


where Is the density of the atmosphere at perigee, h^ Is 
the height of perigee, H Is the height of the homogeneous 
atmosphere. For small orbital eccentricity formula (2) can be 
represented as: 

y ^ c<rp[>/(i-wiu>. tt 


Here: 

p^ Is the density of the atmosphere at apogee, t^ Is the time 
It takes for the satellite to pass through perigee. 

Let us consider equations (1) when: 
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U-;?. . . 


these equations can be written as: 


vj, -(/ t ' A/'W 0,0, 

» -■ ■< o j siVto- > - t^^coiS, 

•i SI 7 > t(fl\(0,CO!,* > 9 ^'.(>t<x)- • 


(3) 


Here the dot indicates a differentiation with respect to t, 


tquacions (3) are simplifications. For example, they do /7 

not allow for so important a fact as the difference in atmos- 
pheric density values at the sunlit orbital segment and that 
in the Earth's shadow. Nevertheless these equations can reveal 
the main features of the dynamics of uniaxial aerodynamic atti- 
tude control. 

We shall point out several properties of equations (3). 

To shorten the notation we shall use the vector conventions. 

We introduce the vector: 


and define the function: 


F(j,t.l^)€ fi* 


6 



ORtGINAL PAGE 19 
OF POOR QUALITY 


such that set (3) can be written as: 

i = Fli.t.if) . (3* ) 

Set (3*) has the property (E) t5] with respect to the 
matrices : 


l.e. the function F(z,T,n) satisfies relations: 

Moreover : 

F(i.r*iyt.f)^F(i.r.f). F(M*^e4.T.^)~S’P{SM,T.f), 

ID/ 

where : 

•( 0 . 0 , 0 . 1 . 0 . 0 )^. S 

When n = 0, set (3') is autonomous: *F(»,t,o)/9t no. 

Set (3) admits of stationary solutions: 

5) , •O', Jlj =a>$8* ; O j ■ - iin Sd ^ •fi ~ O; ■y.st, ^ • (6) 

that describe the position of equilibrium of the satellite in 
the orbital coordinate system. In all these positions of 
equilibrium the itxls coincides with the axis 0 X 2 (Ox^) 

coincides with the axis OX 2 or OX^. 

3. The Integral Surface 

We shall call the movement of the satellite in which the 
angle 6 between axes Ox^^ and (6 • arccos (cosacos3) ) does 
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not exceed a certain value L the regime of uniaxial aerodynamic 
attitude control. For oxanqple, we may take L ■ 15*. Depending 
on the values of the satellite parameters, different approaches 
are possibl« realization of the regime of aerodynamic atti- 
tude ccatrol and its investigation. In particular, we may con- /8 
aider one of the equilibrium positions (6) to be the nominal 
undisturbed motion of the satellite in this regime. In this 
case analysis of the regime of aerodynamic attitude control will 
primarily consist in investigating the stability of the chosen 
position of equilibrium. Such an approach is natural when the 
aerodynamic and gravitational moments acting on the satellite 
are conparable in magnitude < ) • But if the aerodynamic 
moment preponderates it is advisable to conduct the 

investigation of the regime of aerodynamic attitude control in 
different fashion. We discuss below one of the possible methods 
of such an investigation. The analysis will be based on the 
infinitesimal t ->/><' . 

We shall write set (3) as: 

fi, ht a** It), 

£0, (.(1 in (7) 

■ i ' * tf/H OffOUl . 

When e ■■ 0 this set has solutions for which: 

.0, - «»$ir , if ^ . - itjfX . ^ . 0^ (8) 

while the variables y and are determined by equations: 

y - J2, . i ), . - (9) 
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Fig. 3 shows the phase picture of the set (9) In the case of 
p > 0.2. The stationary solutions (6) of set (7) at e ■ 0 
correspond to the stationary solutions of this sett 

Ji-o, 


Motions of type (8) and (9) can be used as nominal undis- 
turbed motions of the satellite In the regime of uniaxial 
aerodynamic attitude control when the aerodynamic moment acting 
on the satellite Is Infinitely large. In these motions t 

.1 fti, 

l.e. the axis Ox^ Is directed exactly along the orbital tangent. 

When e 0, set (7) does not have solutions (8) , (9), but when 
e << 1 It is possible to construct Its Integral surface » similar ^ 
to the family of such solutions, in the form of formal power 
series. The motions that belong to this Integral surface shall 
be regarded as nominal undisturbed movements of the satellite 
In the regime of uniaxial aerodynamic attitude control. 

We shall seek the Integral surface of the set (7) that 
merges Into the family of solutions (8) , (9) when e > 0 In the 
form: 

'' . " * f ) -;ui .1 o,.r. ^ 

/ /**('. / . {/. t ) y . 

where : 

Q, (11) 

We shall regard the above series as formal, l.e. we are not con- 
cerned about their convergence. 
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Inserting ssrist (10), (11) into set (7) end equating the 
expressions for identical powers of e on the left and right, we 
obtain a sequence of recurrent relations: 


A*" * ^ • 

s> u . V * sr ' J'** • 

Q» “ y* ( W j* <•< s) sitt > > • y*4 

(k ... ». 


Here g^j^ (j ■ 1 


,6) are Icnown functions that depend ont 


' • ' . ‘•’f . 

,,.»o ./.V 


when ^ • * ' 

It is easy to see that this sequence uniquely defines the 
un)cnown coefficients of series (10), (11). The coefficients 
of c in these series have the fomt 


s , 



V 


f 

’ ’ q 

.<’s3 

kill H ^ ///, f \/,/ 4/v A t ). 

V 

7 ■ 

• 


7 


\t0tM 



/'^ . ^ ,shii . 


Let us point out several properties of the integral sur- 
face (10) , (11) . To shorten the notation we introduce the 
vector function: 
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% 


y, y(». Sit, ^ >• ^ 


Using relations (4) we can prove that: 

*. • .'/./) - f(S,S),, f ). 

f'( >. v,. i .n = P(s.o,. f.t/.n. 
Q< >. < ) - - •• ) ; 


(12) 


i s 'i'l V.- S', . ■' . 7 . f ) yV'V. ' y. f 

/ ' V. V-.y. c) ■ - /VA, ^). 

(J(x. -V,. 0(X, Vi,T,^>.C). 


(13) 


By virtue of relations (5) we have: 

f(S, V,,r ' ,'v . V- ^ 

V. ’ 

£J(^. li,. ) -Ol'jr. (14) 

y (* f , r. -5^. £ ) - S V(f^. w. ,T. ^ e* . 

I'i ^ . Q,.'f j/.f ) /Y'L (15) 

(1^ ' V/ . 7,y. t ) ^ Vf.ny.f ). 

We can show that when K (>■ ^/z , : 

,f • 'f'. ^ *■*■’* ^<C' • t ^ ) » 

PlK».n.r.i^.n OdcO.T.^.O^o. 


i.e. the stationary solutions of (6) lie on the constructed 
integral surface. 


When n = C. series (10), (11) do not contain x, and the 
integral surface specified by these series can be numerically 


11 


ORIGINAL PAGE ?S 
OF POOR QUALITY 


investigated rather easily. Let us consider the case n ■ 0 in 
more detail. With precision down to terms of order 0 (e) we can 
write set (11) in form (9). This set is equivalent to the 
equation of mathematical pendulum: 

and is integrated in elliptical functions. Let us consider the 
periodic oscillatory and rotational^ solutions of set (11) 

(n = 0) » similar to the periodic oscillatory and rotation solu- 
tions of set (9). Without restricting the generality we shall 
consider that u '' 0 and y(0) “ 0* 

Let us first deal with the oscillatory periodic solutions. 
It is easy to prove that the T-periodic solution of set (9) 
>(t)» for which >(0) * 0, satisfies relations: 

X T ) X'? ' J) -iH). 

0,( 0 0,(r). U,(t ♦ J)- (16) 


by virtue of the last two equations, i2j^(T/4) =0. For the set 
(11) (ti = 0) let us examine the boundary- value problem: 

(To) (17) 

Using the properties (12) , (13) of functions P and Q, we can /12 

prove 15,61 that any given solution >(i), iJj^(t) of this problem 
satisfies relations (16) and, consequently, is T-periodic. For 
the corresponding solution of set (3*); 


^In the case of ti «= 0 we shall call the solution 

of the set (11) ((3*)) the rotational periodic solution if there 

exists a number J'lo (period) such that :(»( t# 

The T values of different sign correspond to rotations in dif- 
ferent directions. 
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(1-0) id) =f |jr(T).o,/r).T. o.« 1 


the equations are valid (cf. (12), (13)): 

Si(-T)« Id). ad ' -J). 


( 18 ) 


by virtue of which: 

i(r^ |)*SS'*d) • 


It follows from this that the second and fifth components of the 
vector z(t) (corresponding to variables ^2 
T/2-periodic functions of t , while the remaining components 
are T/2-antiperiodic f\inctions* Setting t = 0 in (18 ) , we 
find the boundary conditions: 

S HO) ‘210), (19) 

which are satisfied by the solution of z(t). The scalar form 
of these conditions is: 

n , d) . i<>) ^ 

It can be shown [5,6] that in the case n * 0 for any given 
solution of z(t) of the boundary-value problem (3'), (19) 
relations (18) are valid and, consequently, this solution is 
T-per iodic. 

Let us turn to the rotational periodic solutions. Any 
given rotational T-periodic solution of set (9) y(x), ^^^(t), 
for which y(0) * 0, satisfies the relations: 

-S(T), 


By virtue of the first two of these relations we have 
y(T/4) = tt/ 2. For the set (11) (n ® 0) let us consider the 
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m ( 21 ) 

It can be proved that any given solution y(t), of this 

problem satisfies conditions (20) and thus is T-periodic. For 
the corresponding solution of set (3*) (n = 0) : 

%[X) - ] 

the equations are valid (cf. (12), (15)): 

Si(-T)^z(T), 2(t+ ry4>. (22) 

by virtue of which: 

i 

Setting t = 0 we find the boundary conditions: /13 | 

i 

I 

(23) j 

which are satisfied by the solution z(t). The scalar forms 
of these conditions is: 

It can be proved that, when n ^ 0, amy given solution z(t) of 
the boundary problem (3*), (23) satisfies relations (22) and, 
consequently, is a rotational T-periodic solution. 

Numerically solving the boundary problems (19) , (23) for 
the set (3*) in the case n = 0, we can construct periodic solu- 
tions in explicit form that lie on the integral surface (10) , 

(11) . If there exists a family of such solutions for T values 
that belong to a certain interval T^^ < T < T 2 , this family will 
form a subset of the investigated integral surface that wholly 
consists of periodic solutions. The subset has dimensions of 
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2 , and its parameters may be the period and phase (in the solu* 
tions of boundary problems (19) and (23) the phase is fixed, 
but due to the autonomy of set (3') when n * 0 this parameter 
can be chosen at will) . 

4. Numerical Investigation of the Integral Surface 
Tlgy ; (11) when n = 0 

A numerical construction of the solutions of the boundary 
problem (3*), (19) in the case n = 0 was done as follows. A 
certain value of the period T was assigned and for this value 
the problem (19) was solved by the zeroing method. As the 
first approximation of the unknown initial conditions 
f2^(0), ^^ 2 ( 0 ), $(0) we used the values: 

where k is the root of the equation: 


(24) 

K(k) is a complete elliptical integral of the first kind. Such 
values of , ^ 2 ( 0 ), B(0) are obtained by converting from the 

initial conditions of the solution of boundary problem (9) , 

(17) to variables ^2 ^ formulas (8). We observe 

that equation (24) has real roots only when /14 

For X = 0.25, M = 0.2, k = 10.30, Fig. 4 and 5 show the 
dependence on T of the initial conditions flj^(O), (^ 2 ^^)' ^(0) 
of the solutions of boundary -value problem (19) (for k 30 
Fig. 4 shows only Qj^(O) as a function of T) . By virtue of the 
second relation (4) , the curves obtained from the curves in 
Fig. 4 and 5 by the transformation: 
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will also specify the initial conditions for the solutions of 
boundary problem (19) . The curves shown in Fig, 4 and 5 con- 
sist of separate pieces. In the scale of the figures the 
breaks between certain pieces are unremarkable and for clarity 
we show them by circles. The Jacobian: 

j ( 25 ) 

dla,(o). sit(v),fiio)] * 

calculated for solutions belonging to an identical piece, has 
an identical sign. The sign of J changes in moving to an 
adjacent piece. 

Fig. 7 and 8 show graphs of the functions: 

)((r),U(T).fi(r) 10^1 *T) 

for several of the calculated T-periodic solutions of set (3) . 

An analysis of these figures reveals that the breaJcs of the 
curves in Fig. 4 and 5 are due to resonances between the slow 
(with frequency '\<h) and fast (with frequency '^/k) oscillations 
of the satellite. The presence of such resonances suggests a 
divergence of the series (10) and (11). This situation is 
typical for periodic solutions of sets of differential equations 
with a large parameter [7,8]^. 

In order to construct the rotational periodic solutions 
that lie on the integral surface (10), (11) when n = 0, the 
boundary-value problem (3'), (23) was solved. The solution was 
done by the method of zeroing. As a first approximation of the 


■^In 17] an independent set of differential equations is examined 
with small positive parameter y that regularly figures in it, 
and periodic solutions of this set with a period of 'vl//y were 
sought. If, in the set of [7], we replace the independent 
variable x •+• /y and make another series of manipulations, we 
can obtain a set with the large parameter l//y, similar to the 
set (3) . 
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unknown starting conditions ^^^(0), ^^0) the value was 

used: 

where k is the root of the equation: /15 

Such values of obtained by changing from 

the initial conditions of the solution of the bound ary- value 
problem (9), (21) to variables ^2 ^ formulas (8). For 

A = 0.25, \i 0.2, K = 10.30, Fig. 4 and 9 show the dependence on 
T of the initial conditions f2j^(0), ^2^0)' 6(0) of the solutions 
of the boundary problem (23) (for k = 30 Fig. 4 shows only 
fij^(O) as a function of T) . By virtue of the second relation 
of (4) , the curves obtained from those in Fig. 4 and 9 by the 
trans formation : 

- -SJ .(O). 1?, (0)- T-*-r . 

will also specify the initial conditions of the solutions of 
the boundary problem (23). The curves shown in Fig. 4 and 9 
consist of separate pieces. In the scale of the figures the 
breaks between certain of the pieces are not notable and for 
clarity we indicate them by circles. The reason for the 
occurrence of these breaks is the same as in the case of the 
oscillations: resonances between the slow and fast motions of 

the satellite. The behavior of the sign of the Jacobian: 


^(Q, (^).y^).P(? )i 
* I Q j<>). f'foT 

for the calculated solutions of problem (23) is similar tc that 
described above for the sign of the Jacobian (25). Fig. 11 and 
12 show examples of the rotational periodic solutions of the 
set (3) . 

In order to investigate the stability of the discovered 
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periodic solutions of (3) , the characteristic equation of the 
corresponding set of equations in variations was examined. 

Using the property of symmetry of the set (3) and of the 
investigated periodic solutions, it can be proven that this 
characteristic equation is reciprocal. Furthermore, by virtue 
of the independence of set (3) in the case n = 0, this equation 
has a root equal to 1 with multiplicity of at least 2. In view 
of these remarks, this particular characteristic equation can 
be represented as: 

(p-l)\p^-a,p*lXp^-a^p+l)^0. 

where and a^ are certain coefficients. If a^ and a 2 are real 
and: 

ia,|<? • 

all the roots of equation (26) will lie on the circle jp] =1 
and the necessary conditions of orbital stability of the 
investigated periodic solution will be fulfilled. Otherwise, 
this solution is unstable. Figs. 6 and 10 show graphs for the /16 

dependence on T of the coefficients a^^ and a 2 for several of 
the calculated periodic solutions of (3). As was found for the 
solutions shown in Fig. 4, 5 and 9, the necessary conditions 
of orbital stedaility are fulfilled for all values of T, except 
for the small neighborhoods of the points of discontinuity of 
the graphs of the initial conditions and narrow zones (with .i 
width AT ^ 0.05) of parametric resonance. The latter are pro- 
duced in the neighborhood of the points specified by the rela- 
tions ^ -2, a 2 -2, a^^ ^ a 2 (cf. Fig. 6,10).^ Since the 

^By an appropriate replacement of the variables, equations (3) 
can be reduced to the Hamiltonian form: the external moment 

acting on the satellite is susceptible of a force function: 

U - 

Therefore, in keeping with the theorem of Krein-Helfand-Lidskiy 
[9] , zones of parametric resonance will occur in the neighbor- 
hood of only some of the points a^^ ^ a 2 » 
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zones of par5unetrlc resonance on the T axis are very narrow, 
it Is difficult to determine them with reliability. Two such 
zones have been found with certainty. In Fig. 6 and 10 the^ 
are shown by the letter n. For the oscillations (Fig. 6) when 
T = 16.16 we have &2 ~ “2.023; for the rotations (Fig. 10) when 
T = 13.86 we have: 

im a, * - Im a, • o. toi. 

Using the first terms of the series (10) for a and 3 we 
can establish that, when n ^ 0, the maximum angle 6^ between 
the axes Ox^ and OX^ in the case of solutions that belong to 
the integral surface (10) , (11) will not exceed: 

As shown by computations, this estimate is rather good for the 
nonresonant periodic solution. According to such estimate, when 
X = 0.25, \i - 0.2, K > 10 and < 2*, we have 6^ < 4*. 

The obtained result suggests the possibility of using nonresonant 
periodic solutions as the nominal undisturbed movements of the 
satellite in the regime of uniaxial aerodyneunic attitude control. 

5. Numerical Investigation of the Integral Surface 

do), (11) l^en TFO ' 

When n 0 the solutions that belong to the integral sur- 
face (10), (11), in general, are not periodic, and the numerical 
Investigation of this surface is complicated. To checlc for the 
existence of such a surface the set (3) was numerically inte- 

** 

grated on large intervals of time. On the segment 0 < t < 200 it 
solutions of this set were calculated with initial conditions 
(cf. (8)): 


For nonresonant solutions of the boundary problems (19) and 

‘23). ISliUax - l«i(0)l- 

** 

For orbits with an altitude above 250 )cm, this interval com- 
prises more than 6 days. 
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ilt(0)*eos}(p, Q,(o)’-sOtg,. =t(c) mfi(p).o. 


The computations were done for X « 0.25, p « 0.2, k » 10, 
n = -0.1 and n = -0.3. The findings are shown in Fig. 13-24. 
On each figure in the plane: 


crosses indicate the points: 

[i((!tin)(n>oetJJr).Sf(2nn))(n*0.>. . .WO). 


Such figures are usually known as stroboscopic pictures. These 
figures also show for each solution the maximum angle 6^ between 
the axes Ox^ and OX^. This angle is found from formula: 

= mO'X ancroi(toid(T)fO%fi(T)) . 

Oil'- ioon 


A comparison of Figs. 13-24 with Fig. 3 reveals a rather 

good agreement between the obtained stroboscopic pictures of 

the solutions of (3) and the phase pattern of (9). A comparison 

by pairs of Fig. 3 with Fig. 13 and 14, 15 and 16, 17 and 18, 

19 and 20 suggests that, for identical parameters Yq and 

when n = -0.1 the agreement between the stroboscopic pictures 

and the phase curves in Fig. 3 is more accurate than when 

n - -0.3. In all the versions of the analysis it was found 

that 9 < 15*’. This result, as well as the results of $4, 

nt 

testify to the possibility of using the satellite motions that 
belong to the integral surface (10) , (11) as nominal undisturbed 
motions in the regime of uniaxial aerodynamic attitude control. 

The above-considered mathematical model of a satellite is 
quite idealized. This is in keeping with the fact that, at any 
rate, it is described by differential equations that can be 
converted to the Hcuniltonian form. Nevertheless an analysis 
of this model can produce substantive results in the dynamics 
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of aerodynaunic satellite attitude control systems. The methods 
used in the work can be employed to investigate a broad class of 
such systems. 
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Fig. 4. (c « 30. 



Fig. 5. Solutions of the boundary- 
value problem (19), k » lo. 
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Fig. 6. Stability of solutions of 
the boundary problem (19) » »c ■ 10. 
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Fig. 7. Nonresonant solution 
of the boundary problem (19) , 

/ Je4g/S , 




Fig. 8. Resonant solution 
of boundary problem (19) . 

Wf/0, » 

/fw« . 


ORIGINAL PAQr 
OF POOR QUALITY 



Fiq. 9. x’olutions of the boundary- 
value problem (23) , k = 10. 



Fiq. 10. Stability of solutions 
of boundary problem (23), k = 10. 














